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We study the application of the exact renormalisation group to a many-fermion 
system with a short-range attractive force. We introduce a boson field to describe 
pairing effects, and take a simple ansatz for the effective action. We derive a set of 
approximate flow equations for the effective coupling including boson and fermionic 
Q\ fluctuations. The initial conditions are obtained by renormalising the interaction to 

fit the scattering length in vacuum. At some critical value of the running scale, the 
numerical solutions show a phase transition to a gapped phase. Standard results are 
recovered if we omit the boson loops. When boson fluctuations are included, we find 
Q_i' that their contributions are significant only in the small-gap regime. 
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Attractive forces between fermions, even if they are too weak to produce two-body bound 
states, play a crucial role in many areas of many-body physics. In a system of fermions, any 

attraction, no matter how weak, can cause such particles to form correlated "Cooper pairs" 
& [ 

at zero temperature. Examples range from superfluidity in liquid helium-3 to colour super- 
conductivity in dense quark matter. The ground state of the system becomes qualitatively 
different through this pairing, as can be seen from the occurrence of a phase transition at 
some critical temperature. The order parameter for this is the energy gap which appears in 
the fermion spectrum. Within such systems we can identify two extreme limits, depending on 
the interaction strength. The weak-coupling regime, with no two-body bound states, mani- 
fests itself in Bardeen-Cooper-Schrieffer (BCS) superconductivity, while the strong-coupling 
regime, with a deeply bound state which approximates an elementary boson, corresponds 
to Bose-Einstein Condensation (BEC). 

In this letter we study pairing within a framework inspired by modern effective field 



theory (EFT). An EFT description of phenomena is intended to be generic, independent 
of details of the underlying theory, and depending only on the degrees of freedom and 
interactions relevant at the energy scale being considered jjj. An interesting, but rather ad- 
10c approach to combining pairing with EFT ideas is the work of Papenbrock and Bertsch 
2|. A more detailed analysis of the same problem, but cast in a slightly different language, 
can be found in the condensed matter literature |3j, and has a long history J4| (see also 
J5|). There is also some closely related work by Weinberg 6|, concentrating on the effective 
potential in a field-theoretic framework. 

Papenbrock and Bertsch J2J use a standard BCS approach with a contact interaction, 
which they relate to scattering observables by a subtractive renormalisation scheme. Their 
results are limited to positive chemical potential but many of their basic equations are 
identical to those of Marani et al. |3|. In the latter work the chemical potential /j, is not 
tied to the Fermi momentum pp, making it possible to study the crossover to BEC, which 
occurs for negative values of \i. In order to find a more complete way to include many-body 
physics beyond the mean field, within an EFT-inspired approach, we have used a method 
that draws together ideas from Ref . M and Refs. \A M . The approach is based on the use 

nfiij 

of the Exact Renormalisation Group (ERG) [3, |8|,|9(. 

The goal of the ERG approach is to construct the Legendre transform of the effective 
action: T[(f) c ] = W[J] — J ■ (f> c , where W is the usual partition function in the presence 
of an external source J |ll|. The action functional T generates the 1PI Green's functions 
for small fluctuations around the ground state, and it reduces to the effective potential for 
homogeneous systems. Instead of trying to evaluate it directly from W[J], one can use an 
artificial renormalisation group flow, created by introducing an artificial gap in the energy 
spectrum for the fields. This depends on a momentum scale k since we can define the 
effective action by integrating over components of the fields with g > k. The RG trajectory 



then interpolates between the classica 
the full effective action (at k — 0) 



action of the underlying field theory (at large k), and 



la M. This method has been successfully applied to a 
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range of problems, from condensed matter physics 11, 12] to particle physics [13, 14]. 

Here we study a system of fermions interacting through an attractive two-body contact 
potential. We take as our starting point an EFT that describes the s-wave scattering of two 
fermions with a T- matrix determined by the scattering length a®: T = — (47r/M)(— l/a — 
ip)^ 1 . Here a positive scattering length corresponds to a system with a two-body bound 



state (and hence repulsive phase-shifts for low-energy scattering) and a negative scattering 
length to one without a bound state. The binding energy gets deeper as a^ gets smaller, 
while the limit a$ —* ±00 corresponds to a zero-energy bound state. 

Since we are interested in the appearance of a gap in the fermion spectrum, we need 
to parametrise our effective action in a way that can describe the qualitative change in the 
physics when this occurs. A natural way to do this is to introduce a boson field whose vacuum 
expectation value (VEV) describes that gap and so acts as the corresponding order parameter 



. This ERG approach is similar to that for the chiral phase transition in Refs. [13, [14 1 . 
At the start of the RG flow, the boson field is not dynamical and is introduced through 
a Hubbard-Stratonovich transformation of the four-point interaction. As we integrate out 
more and more of the fermion degrees of freedom by running k to lower values, we generate 
dynamical terms in the bosonic effective action. 

For simplicity we treat a single species of fermion, as in neutron matter. We introduce 
a boson field <fi describing correlated pairs of fermions, which leads rather naturally to the 
following Ansatz for V: 



r[ij),ij)\ 0, </>,//, k] = d A x 






(1) 

The field <fi has a non-standard normalisation, since we shall relate its VEV to the gap A 
via A 2 = (0V)- This means that we have only a dimensionless coupling- const ant renormal- 
isation Z g for the boson-fermion coupling. The bosons carry twice the charge of a fermion, 
and so couple to the chemical potential via a term 2ixZ^(f) which has been absorbed into 
the definition of the quadratic term in the potential U. (Here M is the mass of the fermions 
in vacuum and m is naively chosen to have the value 2M, but its real role is only to make 
Z m dimensionless.) 

We expand ^potential U anont its ndninnnn to quart ie order in the Seid flQ, 

U(<f>, 0t) =u + Ul (0V - A 2 ) + ~ « 2 (0V - A 2 ) 2 , (2) 

where the u n are defined as the derivatives of U at its minimum, which occurs where <ft = 
= A. The values of the u n run with the scale of the regulator, k. They also depend on 
the chemical potential fi, the other external parameter in I\ 



In the RG evolution, we start at high k from a free bosonic action and gapless fermions. 
In this symmetric phase we have U\ > and A = 0, and so the VEV of <fr is zero. When k 
is lowered, we expect u\ to decrease until it reaches zero. At this point there is a transition 
to a phase with spontaneously broken U(l) symmetry and a fermion energy gap. Here the 
minimum moves away from A = and, since we expand around the minimum of U, we take 
Mi = in Eq. (J2J). The bosonic excitations in this phase are gapless "Goldstone" bosons. 

The other parameters in the action, namely the wave-function renormalisations Z^, the 
kinetic mass renormalisations Z m ^M and the coupling Z g , are all evaluated in the background 
field corresponding to the minimum of the potential. These other parameters also run with 
k. For large values of k, the fermionic parameters Zm and Z^ tend to unity, while the 
bosonic ones Z$ and Z m tend to zero. We choose the normalisation of the field <fi so that Z g 
tends to unity for large k. The chemical potential /i can be determined from dT/dfi = n, 
where n is the baryon-number density. To be able to study the whole range of phenomena 
from BCS to BEC (where /z becomes negative) we keep this density fixed, and allow fi to run 
during the RG evolution. Following Refs. J3, |^| we define a "Fermi momentum" pp in terms 
of density: Pf = (S^n) 1 / 3 . In the symmetric phase the fermion spectrum does not change 
and pf can be related to [i by fi = p F /(2M). However once a gap appears this connection 
is lost and indeed the whole idea of a Fermi surface may lose its meaning. 

In the results presented here, we allow only Z^ and the parameters in the potential to 
run independently since this is the minimal set needed to treat the bosons dynamically. We 
freeze the other coefficents (Z^, Zm and Z g ) to unity or, in the case of the boson kinetic 
mass, we have also explored setting Z m = Z^. 

The evolution equation for T in the ERG has a straightforward one-loop structure J9(. 
For constant fi it can be written 



d k T = -- Tr 



BB 



I 



Trr — R-b) <9fcR-B + 77 Tr (T FF — R_p) <9fcR 



2 



,(2) 



FF 



(3) 



(2) 

Here ^ F f(bb) ^ s ^ ne ma trfx containing second functional derivatives of the effective action 
with respect to the fermion (boson) fields and Rb(b) is a matrix containing the corresponding 
boson (fermion) regulators. A 2 x 2 matrix structure arises for the bosons because we treat <j) 
and $ as independent fields in order to include the number-violating condensate. A similar 
structure also appears for the fermions. By inserting our ansatz for V into this equation we 
can turn it into a set of coupled equations for the various parameters discussed above. 



The one-loop structure of the ERG equations means that the evolution at large k can be 
matched onto that of the underlying EFT 1], albeit for a non-standard choice of regularisa- 
tion. In the bosonic sector, the regulator is an additional bilinear term in the fields. It has 
the matrix structure 

R B {q,k)=R B {q,k)dmg{l,l), (4) 

where Rs{q,k) is a scalar function. This regulator provides an extra contribution to the 
single-particle energies, which should suppress the contributions of states with momenta 
q ^ k. As k tends to zero Rsik) should vanish, so that we would recover the full effective 
action in the absence of any truncations. For large k, Rs{q, k) should be large for momenta 
q ;$ k, to give all these modes large energies. In particular, it should be of order k 2 if the 
behaviour of any integral for large k is to reflect its order of divergence. For more discussion 
of the choice of regulator in this approach, see Ref . |l5| . 

In the symmetric phase, our regulator for the fermions should be positive for particle 
states (q > p F ) and negative for hole states (q < p F ), so that it provides a /c-dependent 
energy gap for excitations around the Fermi surface. It should have the structure 

R F (q,p F , k) = sgn(e(g) - /i) R F (q,p F , k) diag(l, -1), (5) 

where e(g) = g 2 /(2M). The function R F (q,p F , k) should be peaked about p F to suppress the 
contributions of states with momenta \q — p F \ ^ k, but otherwise it should behave like the 
bosonic function. Note that the "Fermi surface" in the gapped phase no longer lies exactly 
at p F , and there may not even be a well-defined Fermi surface. However, in that case a real 
gap has appeared in the spectrum and the regulator no longer plays a crucial role. 

It is worth noting that we can recover the usual mean-field results if we include fermion 
loops only, ignoring all diagrams with virtual bosons. In this case the equation for the 
effective potential can be integrated analytically without making any truncation, and the 
results for k = coincide with those in Refs. J2|, |3| . The initial conditions on the evolution 
are obtained by assuming that in vacuum our theory reproduces the scattering length a , 
and that the necessary subtraction are identical in matter and in vacuum. (See Ref. yj and 
references therein for discussion of the vacuum problem.) The resulting running potential 
has the form 



where we have introduced the short-hand notation 

E FR (q, k) = e(q) - \x + R F (q,p F , k) sgn(e(g) - //). (7) 

For k = this has a closed-form expression in terms of an associated Legendre function, 
PTiv): 

W,„,o> = J£- (j-L- - -L(, + V)^ (-^)) ■ (8) 



where &a = v2MA and xo = ///A. In the weakly attractive limit, p F ao — > 0~, the chemical 
potential is directly related to the density by p = e F = p%I(2M), and minimising JJ MF with 
respect to A leads to the standard result (see, e.g., Ref. 2j) 



A = -e F exp(--^ r -). (9) 

More generally, working at constant density, we can minimise U subject to the subsidiary 
condition dU /dp = n. This gives the results from Ref. jj|, which can be turned into a pair 
of equations for A and p, 

A/e F = (-37r 2 (^ + l)^pi /2(2/0 ))- 2 / 3 

1 2 7 1 \ / A \ V2 






g + 1) 1/4 (^2 (ito) - i ^/2 (Ito)) (^) (10) 



p F a 3 

where yo = -^0(^0 + 1) -1 ^ 2 - Note that, in this case ,A/e F can be eliminated from these 
results to leave a single equation for xq = ///A in terms of p F a,Q. 

We now turn to the full evolution equations. At the present level of truncation, all of 
these can be obtained from the evolution of the effective potential. It is convenient to write 
this as U(p, /i, k) where p = <p*<p. The density and Z^ are then given by 



dU_ 
dp, 



z l d2u 

p=A2 ' ~" 2 dpdp 



(11) 

p=A 2 



By evaluating the loop integrals in Eq. (J3j) for a uniform boson field, we get evolution 
equation for U at constant p,, 

+ J_fjPgE ss= 



2Z^ J (2vr)3 v ^|^ _ v; 
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where 

E BR {q) = ^-q 2 + Ul + u 2 {2p-A 2 )+R B {q,k), V B = u 2 p, (13) 

and Efr is defined above. Substituting our expansion on the left-hand side leads to a set 
of equations for the u n , each of which is coupled to the coefficient of the next term through 
the running of A 2 , the position of the minimum of U. These equations have the form 



d k u n - u n+1 d k A' 



Qn 

~djf l 



d k U 



(14) 



p=A 2 



In the symmetric phase (A 2 = 0) they decouple. In the broken phase they do not and, for 
example, the left-hand-side of the equation for d k u 2 contains the coefhcent u%. We could 
simply set u^ = 0, but a better approximation can be obtained by substituting the form for 
Us(k) from the evolution with fermion loops only, as described above. This is the approach 
we adopt here for w 3 and similar coefficients in the equations for n and Z^. It has the 
advantage of providing an approximation that becomes exact in situations where boson 
loops can be neglected. 

In order to follow the evolution at constant density we allow \x to run with k and define 
the total derivative 



d k = d k + (d k p) 



d_ 
dp 1 



(15) 



where d k p = dp/dk. Applying this to dU /dp and demanding that n is constant (d k n = 0) 
gives 

-IZ^ 4 A 2 + X d k p = - — ( 0,1 



~dJl( bkL 



=A 2 



where we have introduced the fermion- number susceptibility x = d 2 U /dp 2 \ p= &2 
tions for the coefficients u n in the potential and Zj, are 



d k u + n d k p 
d k ui - u 2 d k A 2 + 2Z (j> d k p 

d k u 2 - u 3 d k A 2 + 2z (f)1 d k p 



dkU\ p=A 2 , 

d k U 
d 2 



U 1 



dkZ^ - z^i d k A + -x'd k p = - 



dp 2 
1 



d k U 
d 2 



where we have defined 



Zt/,1 



1 d 3 U 

2 d 2 pdp 



X 



p=A 2 



2 dpdp 
dp 2 dp 



p=A 2 
p=A : 

d k u 



p=A 2 



(16) 
The equa- 

(17) 
(18) 

(19) 

(20) 

(21) 



p=A 2 



8 

The driving terms in these evolution equations are given by appropriate derivatives of 
Eq. ([12)1 . In the symmetric phase we evaluate these expressions at A 2 = 0. The driv- 
ing term of Eq. (jl6|) vanishes in this case, and hence [i remains constant. In the broken 
phase we keep A 2 non-zero and set u\ = 0. 

On the left-hand sides of these equations, the coefficients M3, z^i, \ an d x' an correspond 
to terms outside our present truncation and so we replace them by their expressions from 
evolution with fermion loops only. This approximation provides a closed set of equations, 
whose solutions agree with the exact results when boson loops are neglected. Underlying it 
is an assumption is that bosonic contributions to the evolution are not too large compared 
to fermionic ones. We shall investigate this assumption below, and find that it seems to hold 
for a wide range of values of pfclq- In future work, we hope to be able to check it further by 
including more bosonic terms in our ansatz for T and examining convergence with respect 
to the number of terms. 

We integrate the resulting differential equations numerically. In this work, we use regu- 
lator functions with the forms 

k 2 k 2 

R F {q,k;p F ,(r) = —6{q-p F ,k;a), R B (q,k;a) = —6(q,k;a). (22) 



where 6 is the smoothed step-function 
9(q,k;a) 



2 erf (I/a) 



erfl-MVerf^ {q ~ k) 



(23) 



ka J \ ka 

with a being a parameter controlling the sharpness of the step. The derivative dk9 falls off 
for k — > 00 and so dkR.B,F provides a UV cut-off on the loop integrals in Eq. (J12J1 [9J. 

We run the evolution down towards k = starting from some large, fixed scale, k = K. 
The initial conditions for this are obtained by matching onto the standard evolution in 
vacuum for k > K using Eq. (J^J) and related expressions. Differentiating U MF with respect 
to A 2 at A = we get the initial condition on Ui, 

„ IK\ M 4. 1 f d ' q f sgn(e(,;) - y.) 1 \ 

The second term in the integral contains exactly the same linearly divergent term as in the 
free inverse T matrix. It cancels with the similar divergence in the first term to leave a finite 
result. This result is linear in K, reflecting the divergence of the underlying integral. It also 
contains a correction for the fact that our regulator depends on the density (through p F ). 
The initial values of u<i and Za, can be obtained similarly from second derivatives of U MF . 



0.15 




k(fm ) 
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FIG. 1: Numerical solutions to the evolution equations for infinite ao and pf = 1.37 fm, starting 
from K = 16 fm -1 . We show the evolution of all relevant parameters for the cases of fermion loops 
only (orange/grey lines), and of bosonic loops with a running Z^ (blue/black lines). All quantities 
are expressed in appropriate powers of fm . 



Although our results can be applied to many systems, for definiteness we concentrate 
initially on parameter values relevant to applications to neutron matter: M = 4.76 fm -1 , 
Pf = 1.37 fm -1 , and large two-body scattering lengths (|a | > 1 fm). We then explore a 
wider range of values for ppcio. We have checked for dependence of our results on the starting 
scale K and find that this is undetectable as long as K > 5 fm -1 (about 4^op). Similarly, we 
get numerically indistinguishable results for a range of values of the width parameter a. 

Some typical solutions for the evolution equations are given in Fig. ^ for the case of 
infinite ao- We compare two different approximation schemes, one where we have fermion 
loops only, and one where we include boson loops as well and we allow Z^ to run. For large 
values of k the system remains in the symmetric phase. At k CT - lt ~ 1.2 fm -1 , u\ vanishes and 
below that point the system is in the broken phase where we plot A instead of u±. 

One immediate observation is that the contributions of boson loops are small. Indeed in 
the symmetric phase their inclusion has essentially no effect. Below the transition they do 
become visible, particularly in U2- However their effects on the gap are even smaller, at most 
~ 1% for l/|jOpao| < 1. In this region, the main effect of the bosons is a small enhancement 
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FIG. 2: Fractional deviations (in %) of the gaps for various numerical solutions from the analytic 
mean- field result. All cases are for pp = 1.37fm~ . Solid dots (red) denote the results from a 
calculation with fermion loops only, the open circles (blue) show the results when boson loops are 
added, and the crosses (green) show the effect of allowing Z^, to run as well. 



of the gap related to the reduction in u 2 as k — > 0. Not plotted are results including boson 
loops but with Z^ fixed at unity; the results are vitually indistinguishable from those with 
running Z^. The same is true if we set Z m = Z^ rather than Z m = 1. 

It is also instructive to examine the behaviour of the gap over a wider range of l/(pj?ao). 
The overall picture the same as found in Ref. [5( for fermion loops only, with a crossover 
from BCS pairing (with positive /i ~ ep) for l/(ppCho) < to BEC (with increasingly 
negative values for /i) for l/(ppao) > 0. For negative 1/(pfOo), the gap lies close to the 
exponential curve of the weak coupling limit, Eq. Q. However a closer look at these results, 
as in Fig. shows that deviations from mean field behaviour are present in this region and 
become increasingly noticable for weaker couplings or lower densities. 

We are not able to follow our results beyond l/(ppa ) ~ —2. This is because of the non- 
analyticity of the effective action for small gaps which means that our expansion Eq. (J2J) is 
no longer adequate. For example, the fermion loops contain a term <j)^<ft\og(<j)^<j)) which gives 
a divergent contribution to «2- Although the complete fermionic part has been incorporated 
in our equations, the analogous boson effects are not included. Some of our results suggest 
that the gap vanishes as a power rather than an exponential, but a definite conclusion will 
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have to await a more complete calculation. 

The dominance of fermion loops for most values of l/(pi?Oo) may seem surprising. How- 
ever this is the region of BEC or the crossover to it, where the mean boson field (the gap) 
is large. In fact the gap is particularly insensitive to boson fluctuations as a result of can- 
cellations between the direct contribution to the running of A 2 and indirect ones via Ui- 
Contributions to other parameters such as u 2 and Z^ are larger, ~ 10% for l/(p F a ) ~ 0. 
Our results imply that these effects are likely to be important for neutron matter, where 
calculations with realistic interactions lead to gaps of at most 5 MeV (~ O.lep) 16}. 

However, if we were to apply our current calculation to neutron matter, we would find a 
gap comparable to e F , of the order of 30 MeV. Fayans has given a simple explanation for 
the smaller values found using more sophisticated treatments [l_3] - 1 The argument can be 
given most succinctly for weak coupling, where the gap satisfies a generalisation of Eq. , 

A = (8/e 2 )e F exp (- (vr/2) cot(6(p F ))) . (25) 

For nucleon-nucleon scattering, cot 5 increases relatively quickly with momentum and the 
resulting reduction in the gap is substantial. We therefore expect that an extension of our 
approach to include the effective range should capture this physics. 

There a number of improvements which could be made to our approach. Adding an 
effective range is clearly an important one. Another is the evolution of the boson kinetic 
mass (Z m ) since this will allow a scaling analysis of boson loops for small gaps. We should 
also include running of the fermion renormalisation constants and "Yukawa" coupling, to 
provide a full treatment of the action Q. Beyond that we would like to treat explicitly the 
particle-hole channels (RPA phonons) since these contain important physics. They will also 
allow us to remove the "Fierz ambiguity" associated with our bosonisation of the underlying 
contact interaction J21[. We plan to explore these extensions in future studies. 
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